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Abstract. The effect of rotation on the cooling of neu- 
tron stars is investigated. The thermal evolution equations 
are solved in two dimensions with full account of general 
relativistic effects. It is found that rotation is particularly 
important in the early epoch when the neutron star's inte- 
rior is not yet isothermal. The polar surface temperature 
is up to 63 % higher than the equatorial temperature. This 
temperature difference might be observable if the thermal 
radiation of a young, rapidly rotating neutron star is de- 
tected. In the intermediate epoch (10^ 10^ yr), when 
the interior becomes isothermal, the polar temperature is 
still up to 31 % higher than the equatorial temperature. 
Afterwards photon surface radiation dominates the cool- 
ing, and the surface becomes isothermal on a timescale of 
~ 10^ yr. Furthermore, the transition between the early 
and the intermediate epochs is delayed by several hundred 
years. An additional effect of rotation is the reduction of 
the neutrino luminosity due to the reduction of the central 
density with respect to the non-rotating case. 
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1. Introduction 

In the past three decades the cooling history of neu- 
tron stars was investigated by several authors (e.g. 
Tsuruta 1966, Richardson et al. 1982, Van Riper 1991, 
Schaab et al. 1996, Page 1997). Recent numerical simula- 
tions account for non-isothermal interior, as well as for 
general relativistic effects. Nevertheless, as far as we know, 
all investigations assumed spherical symmetry of geome- 
try and temperature distribution. 

As it was pointed out by Miralles et al. (1993), the 
effect of rapid rotation on the cooling of neutron stars 
can be as important as general relativistic effects, whereas 
the effect of slow rotation should be negligible. Although 
the assumption of slow rotation holds for most of the 
known pulsars, there exist a couple of millisecond pul- 
sars (Lyne 1996), for which rotation should yield a rather 
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different cooling behaviour. These millisecond pulsars are 
generally located in binary systems. It may however be 
expected that young, isolated millisecond pulsars can be 
detected in the near future, too. A candidate might be the 
supernova remnant of SN1987A. Although the observed 
neutrino burst lasting for ten seconds indicates that a 
neutron star was formed in the supernova, there is no ev- 
idence for the continued existence of it (s. Chevalier 1997 
for a recent review). However, the neutron star can still 
be hidden by the surrounding matter, and the continued 
observations might reveal a rapidly rotating neutron star. 

The aim of this letter is to study the effect of 
non-spherical geometry on the cooling of neutron stars. 
As far as we know, this is the first investigation of 
rotational effects beyond the isothermal core ansatz 
(Miralles et al. 1993) and also the first completely two di- 
mensional simulation of neutron star cooling. The letter 
is organized as follows: We first derive the general rela- 
tivistic equations of thermal evolution and describe the 
numerical method in Sect. 2. In Sect. 3, we apply the two 
dimensional cooling code to static and rotating neutron 
star models based on a relativistic equation of state in- 
cluding hyperonic degrees of freedom. Finally, we summa- 
rize our conclusions and discuss further improvements and 
applications of the current work in Sect. 4. 

2. Equations and Numerical Method 

Already a few seconds after the formation of a neutron 
star in a supernova, its interior settles down into catal- 
ysed, degenerate matter. The subsequent cooling involves 
only thermal processes and does not change the space 
time geometry. However, the structure of a neutron star 
depends on the rotational velocity, which generally de- 
creases as the star loose angular momentum, e.g. due to 
emission of magnetic dipole radiation. Since the time scale 
for reaching hydrostatic equilibrium is much smaller than 
the time scale for the variation of angular velocity (s. 
Zeldovich & Novikov 1971, p. 239), one can treat the evo- 
lution of a neutron star in quasi stationary approximation. 
Though the partial transformation of rotational energy 
into thermal energy may considerably change the cooling 
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behaviour of a neutron star (see, e.g.. Van Riper 1991), 
and also the variation of space time geometry might have 
an eflFect on it, we study here, as a first step, the simplest 
case of constant angular velocity. 

The stationary, axisymmetric, and asymptotic flat 
metric in quasi isotropic coordinates reads 



ds^ 



e^^dt^ + e">'{dip - N-^dty + e^^idr^ + r^d0^) 



(1) 



where the metric coefficients = g^^ (r, 0) are functions 
of r and 9 only. The metric coefficients are determined by 
the Einstein equation (c = G = 1) G = SttT, and the 
energy-momentum conservation V • T = 0. The obtained 
elliptic differential equations (Bonazzola et al. 1993) are 
solved via a finite difference scheme (Schaab 1998) once, 
before the cooling simulation starts. 

In the case of uniform rotation, = const., con- 
sidered here, the equations for thermal evolution are 
(Miralles et al. 1993, Schaab 1998) 



drhi + ^deh = -re'^+2'^ (^ie'^^e + TCydtT 

drf = hi 
rK 

-def = -—e-'-1'h2, 
r rK 



where 



r= (l -e^f'^-''' {Q-N-^fy 



T = -eT. 



(2) 

(3) 
(4) 

(5) 

(6) 

(7) 



h denote the heat flux 3-vector in the comoving frame, 
Cy the heat capacity, e the neutrino emissivity, and k the 
heat conductivity. The partial radial and angular differen- 
tials are abbreviated by dr and dg, respectively. Thermal 
equilibrium is described by T = const.. 

At the surface of the neutron star the heat flux hi and 
/i2 is determined by the normal heat flux /in 



(8) 



/i2(r = R) = hl 



1 fdRY\ 1 dR, 
^+R^[J9 Ed^^^' 



(9) 



where R{9) is the r-coordinate of the surface, /in is 
taken from a non-magnetic photosphere model which 
describes the temperature gradient in the region be- 
tween e = 10^°gcm~^ and the star's surface (e.g. 
Gudmundson et al. 1983). In these models /in(^) depends 



on the temperature at the density e = 10^" gem and on 
the surface gravity 



1/2 



surface 



(10) 



The parabolic differential equations obtained after in- 
serting Eqs. (3) and (4) into Eq. (2) are solved via an 
implicit finite difference scheme by using a alternating di- 
rection implicit method. This yields a non-linear equation 
system which can be solved iteratively. The obtained lin- 
ear equation systems have tridiagonal coefficient matrices 
which can be inverted rather fast. The correctness of the 
two dimensional code was checked by comparing the out- 
come of it with simple, analytically solvable models and 
with the results of the one dimensional code described by 
Schaab et al. (1996). 

3. Results 

We consider a superfluid neutron star model basing on the 
relativistic Hartrce-Fock equation of state labelled RHF8 
in Huber et al. (1997), which accounts for hypcronic de- 
grees of freedom. The global properties of uniformly ro- 
tating models with fixed gravitational mass M = 1.5 Af© 
and angular velocity il = 0, 0.5, and 0.99 i^K are summa- 
rized in Tab. 1. fix denotes the maximum possible Kepler 
angular velocity, above which mass shedding sets in. All 
models allow for both the direct nucleon Urea and for the 
direct hyperon Urea processes (cf. Prakash et al. 1992). 
All direct Urea processes are suppressed by nucleon 
and lambda pairing below the respective critical tem- 
perature (cf. Schaab et al. 1998). The ingredients to 
the cooling simulations are similar to those discussed 
by Schaab et al. (1996) in detail and are published 
on the Web (http://www.physik.uni-muenchen.de/sek- 
tion/suessmann / astro / cool/schaab.0 1 98 /input .html) . 

Figure 1 shows the evolution of the surface tempera- 
ture as measured by a distant observer. It is possible to 
distinguish three epochs of evolution. In the first epoch t < 
100 yr, large temperature gradients occur in the interior 
of the neutron star. Radial temperature gradients were 
already found in one dimensional simulations (see, for ex- 
ample, Richardson et al. 1982). In rotating neutron stars, 
azimuthal temperature gradients, which cause transverse 
heat fiow, /i2 7^ 0, exist, too. The polar temperature is by 
16 % or 63 % higher than the equatorial temperature for 
the models rotating with = 0.5 Qk and 0.99 J^k, respec- 
tively (see also Fig. 2). After about 100 yr the cooling wave 
reaches the surface and the interior becomes isothermal, 
T = const. The temperature deviation between pole and 
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Table 1. Models of uniformly rotating neutron stars with fixed 
gravitational mass M = l.SM© and angular velocity f2 = 0, 
0.5, and 0.99 Ok 



- Q=0.5 

- fl=0.99 



Model 


n = 


n = 0.5 fix 


n = 0.99 Qk 


M [Mq] 


1.5 


1.5 


1.5 


Mb [Mq] 


1.694 


1.693 


1.672 


RT^ [km] 


12.81 


13.25 


17.03 


n [s-i] 





2631 


5208 


ec [lO^^' gcm-^] 


11.67 


10.84 


7.208 


ric [fm"'^] 


0.6285 


0.5888 


0.4065 


g^'^ [lO^cm^s-i] 


1.648 


1.459 


0.546 


51" [lO^cm^s-i] 


1.648 


1.630 


1.396 



Entries are the gravitational mass M, the baryonic mass Mb, 
the circumferential radius as measured by a distant observer 
R^, the angular velocity Q, the central mass density ec, the 
central baryon density ric, and the surface gravity Qs- 




Fig. 1. Thermal evolution of rotating, superfluid neutron stars. 

The thick (thin) curves correspond to the equatorial (polar) 
surface temperature as measured by a distant observer. The 
considered models are given in Tab. 1. 



equator is now mainly caused by the smaller surface grav- 
ity gs at the equator, since Tg oc for fixed internal 
temperature. Whereas the temperature deviation is still 
high in the case of nearly Kepler rotation 31 %), it al- 
most disappears for the model with — 0.5 fix 4%). 
After about 10^ yr, the photon surface radiation domi- 
nates the cooling. The temperature distribution tends to 
a new equilibrium with dTiii((?) /dt = const., on a timescale 
of ~ 10^ yr. T^{9) denotes the temperature at the inner 
boundary of the photosphere at e = 10^°gcm~^. Since 
the heat capacity is only weakly temperature dependent 
for T < 10^ K in the outer crust, the surface temperature 
tends to an isothermal state, whereas the temperature in 
the outer crust varies with the surface gravity. This be- 
haviour is accompanied by the breakdown of the scaling 
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Fig. 2. Fraction /T^'^ of polar and equatorial surface tem- 
perature as measured by a distant observer for the two rotat- 
ing models. The corresponding relations {gs'^/gf)^^^ of surface 
gravity are shown by the horizontal lines. The thin dashed 



curve is calculated under the assumption that Tg oc g, 
even for small surface temperatures. 



1/4 



holds 



Tg OC (^y^ of the surface temperature with surface gravity 
for small surface temperatures. However, as one can see by 
comparing the thin dashcid curve, which assumes Tg oc 
over the whole temperature range, with the thick curve in 
Fig. 2, this is only a small effect and cannot explain the 
tendency of the surface temperature fraction to unity. 

Besides these effects on the angular dependency of the 
temperature, rotation has also a net effect on cooling in 
the intermediate epoch. This effect is caused by the reduc- 
tion of the neutrino luminosity, which sensitively depends 
on the central density. A further effect of rotation on cool- 
ing of neutron stars is the lengthening of the thermal diffu- 
sion time by several hundred years. Additionally the drop 
of the surface temperature is smoother in the case of rota- 
tion (see also Miralles et al. 1993), since the cooling wave 
reaches the polar region earlier than the equator. During 
the cooling wave is reaching the different surface regions, 
the fraction TP/T^'i falls below unity (s. Fig. 2). 

4. Conclusions and Discussion 

In this letter, we have studied the effect of non-spherical 
geometry on the cooling of neutron stars. Our general find- 
ing is that rapid rotation plays a significant role in the 
thermal evolution. This role is particularly important in 
the early epoch, when the star's interior is not yet isother- 
mal. It is precisely this epoch, when even isolated pul- 
sars might rotate very rapidly. The angular dependent ra- 
dial and transverse heat flow cause azimuthal temperature 
gradients in the interior of the star and thus also at the 
surface. The polar surface temperature is by up to 63% 
higher than the equatorial temperature, The simulation 
of this non-isothermal epoch cannot be performed with 
isothermal core approaches or one dimensional codes. 

In the intermediate epoch, heat conduction is unimpor- 
tant, since the interior is nearly isothermal. Nevertheless 
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two significant, though smaller, effects remain. The first 
effect is the reduction of the central density with respect to 
the static, non-rotating case, if one fixes the gravitational 
mass or the rest mass. Since the neutrino emissivity de- 
pends sensitively on the density, especially in the vicinity 
of threshold densities for fast neutrino emission processes, 
the neutrino luminosity of the star is reduced as well. The 
second effect is the angular dependency of the surface tem- 
perature caused by the angular dependent surface gravity. 
The polar surface temperature is still by up to 31 % higher 
than the equatorial temperature. Both effects could also 
be studied with one dimensional codes or even within the 
isothermal core ansatz (see Miralles et al. 1993). Whereas 
Miralles et al. (1993) assumed that the interior is still 
isothermal in the late photon cooling epoch, we find that 
the surface temperature distribution itself tends to an 
isothermal equilibrium, which implies that the tempera- 
ture distribution in the outer crust becomes angular de- 
pendent. 

The transition between the early and the intermedi- 
ate epochs turns out to be delayed by rotation. Because 
of the enlargement of the crust in the equatorial plane, 
the cooling wave formed in the inner region of the star 
needs several hundred years longer to reach the surface. 
Furthermore, the transition becomes smoother, since the 
radial diffusion time it not isotropic any more (see also 
Miralles et al. 1993). 

The surface temperature, especially in the intermedi- 
ate epoch, is rather sensitive to micro physical parame- 
ters (superfluid energy gap, neutrino emissivity, equation 
of state, etc.; see Schaab et al. 1996 for a recent review). 
The comparison of the net effect with observation might 
therefore be difficult. Nevertheless, the comparison of the 
observed soft X-ray spectra with the theoretical obtained 
spectra may reveal the angular dependency of the surface 
temperature due to rotational effects. The deduction of 
the spectrum from the surface temperature distribution 
is, unfortunately, not at all trivial in non-spherical geom- 
etry (Cunningham 1975). So far, only the thermal spectra 
of slowly rotating neutron stars (fi = 70.37s~^ for the Vela 
pulsar 0833-45) are known. It may however be expected 
that the thermal spectra of a young, isolated millisecond 
pulsar can be detected in the near future. For example, 
the period P = 1.56 ms of the fastest pulsar B1937-I-21 
(Backer et al. 1982) known today corresponds to a ratio 
= 0.77, which lies between the respective ratios of 
our two rotating models. 

Since we neglected some important effects, like differ- 
ential rotation and the loss of angular momentum, our 
work can only be understood as a first step towards two di- 
mensional simulation of cooling of rotating neutron stars. 
Of course, the two dimensional code presented here may 
also serve to investigate other non-spherically symmetric 
effects, such as the effect of strong magnetic fields on heat 
conductivity and neutrino emissivity. The inclusion of dif- 
ferential rotation and strong magnetic fields, as well as 



the mentioned deduction of the photon spectra will be 
addressed to future work. 
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